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£h . Abstract 



A correction to the nuclear functional is proposed in order to im- 
prove the density of states around the Fermi surface. The induced 
effect of this correction is to produce a surface-peaked effective mass, 
[~^. | whose mean value can be tuned to get closer to 1 for the states close 

to the Fermi energy. In this work we study the effect of the correction 
^"•O | term on global properties of nuclei such as the density of states in 40 Ca 

and 208 Pb, pairing and specific heat at low temperature in 120 Sn. In 
the latter application, an explicit temperature-dependent form of the 
correction term is employed and it is shown that the critical tempera- 
ture is reduced by 40-60 keV. 



In finite nuclei, the single particle states around the Fermi energy are 
known to be strongly affected by the dynamical particle- hole correlations [1] . 
The energy of the single particle states is modified and thus the level den- 
sity around the Fermi energy is changed, which in turn has an impact on 
low-energy properties such as pairing correlations collective modes [4] 

as well as on temperature-related properties such as the entropy, the critical 
temperature, the specific heat [5]. In nuclear astrophysics, a good descrip- 
tion of the density of states around the Fermi energy (which is related to 
the nucleon effective mass) turned out to be relevant also for the energetics 
of core-collapse supernovae [HE]. 
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Despite the important role of the density of states in self-consistent mean 
field theories, most of these models such as those based on Skyrme [8], 
Gogny [9] and M3Y interactions [10J or on the relativistic approaches like 
RMF [X 1 J or RHF [12] . have a density of states around the Fermi energy 
that is too low. Already in the 1960s, G.E. Brown et al. suggested that 
the effective mass m*/m should be close to 1 around the Fermi energy 
in order to improve the level density [13] . It could indeed be shown that 
the first order expansion of the energy-dependent self-energy around the 
Fermi surface (bHfT5] produces an effective mass which is the product of two 
different terms, the &-mass and the w-mass m w [1]. The w-mass is related 
to the energy dependent part of the self-energy, dynamically generated by 
the coupling of the particles to the core-vibrations |16H19| . These dynamical 
correlations, which lead to an increase of the effective mass at the surface, 
are implemented beyond the mean field. Indeed, the effective mean field 
theories have an average effective mass m*/m, the A;-mass, around 0.6-0.8. 

The different microscopic calculations of the particle-vibration coupling 
(PVC) have been performed at the first order in perturbation, due to their 
heavy computational features (T6H22] , However, the induced effects of the 
PVC on the mean field itself is non-negligible. Since the effective mass is 
enhanced at the surface of nuclei [20] , it impacts the single particle states 
and the pairing correlations which, in turn, modify the PVC and the effective 
mass. There is then a self-consistent relation between the properties of the 
single particle states around the Fermi energy and the PVC. 

In this paper, we propose an effective short-cut for treating the effects 
of the PVC directly in the Energy Density Functional (EDF) approach . 
The coupling of the collective modes to the single particle motion induces 
a dynamical type of correlation that in principle could not be easily imple- 
mented in an effective nuclear interaction or in a nuclear EDF. However, the 
core-vibration is mainly located at the surface of the nuclei which makes the 
implementation in energy-density functionals easier. A parametrization of 
the w-mass as a gradient of the nuclear profile has shown to give good results 
within a nuclear shell model [15J. In addition, the first order expansion of 
the self-energy near the Fermi energy induces a renormalization of the single 
particle Green function as well as a correction to the mean field, which al- 
most compensate the effective component of the equivalent potential |15| . In 
the present approach, the surface-peaked effective mass is included through 
a correction term in the Skyrme energy density functional. This correction 
is energy independent and designed in such a way to have a moderate effect 
on the mean field. 

The paper is organized as follows: in Sec. [1] we will describe the theoret- 
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ical framework. In Sec. [2] we will discuss the application to spherical nuclei, 
in particular, to 40 Ca and 208 Pb. In Sec. [3] we discuss the superfluid prop- 
erties in connection with the surface-peaked effective mass, both at T = 
and at finite temperature. Finally, in SecSJwe will give our conclusions and 
outlooks. 



1 Nuclear energy density functional 

The energy density functional derived from the Skyrme interaction will be 
considered as follows. The standard Skyrme energy-density %(r) is ex- 
pressed in terms of a kinetic term /C(r) and an interaction term written 
as [23] 

H(r) = £(r) + Y, n T(r) , (1) 

T=0,1 



with: 



h 2 

/C(r) = — r(r), 



2m 

H T (r) = C p T p 2 T {v) + Cf V(r)V 2 P T(r) + C t tPt (v)t t (v) 

+ C^ 2 T {v)+C^ J p T {v)V ■ J T (r) , (3) 

where the kinetic density r(r), the density /?r(r), the spin-current JT^(r) and 
the current J(r) as well as the relation between the Skyrme parameters and 
the coefficients in Eq. ([3|) are defined in Ref. [23j. We have kept only the 
time-even component for the application considered in this paper. Notice 
the difference between the integrated energy, H, and the energy density 
functional, 7i. These two quantities are indeed related as: H = J dr7i(r). 

In the following, we explore the impact of adding to the functional an 
iso-scalar correction of the form: 

?C rr (r) = C T(Vp)2 r(r) (Vp(r)) 2 + C^V) 2 (Vp(r)) 2 , (4) 

where the first term have been first introduced in Ref. [28P29] and induces a 
surface-peaked effective mass while the second term is introduced to mod- 
erate the effect of the first one in the mean field. 

The effective mass is obtained from the functional derivative of the en- 
ergy H and is expressed as (q runs over neutrons and protons: q = n,p): 

^ = ir = f + + ^ (Vp)2 (Vp(r)) 2 , (5) 

2m* (r) 0T q 2m H 
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where the coefficient C q = (Cg ± CD/2 (with + for q = n and — for q = p), 
and the mean field reads: 

U q (r) = U^(r) + U c ^(r) , (6) 

where J7 g ky (r) is the mean field deduced from the Skyrme interaction |23j 
and C/ corr (r) is the correction term induced by Eq. dH) which is defined as: 



U co "(r) 



-2C, 



.r(Vp) 5 



T(r)V 2 />(r) + Vr(r)Vp(r) 



- 2C, 



P 2 (V P ) 5 




p(r)(Vp(r)f + /9 (r) 2 V 2 p(r)) . 



(7) 



From the variation of the total energy H with respect to the ground 
state density matrix we obtain the following set of self-consistent Kohn- 
Sham equations: 

h 2 

-V + U q (r) -iW q (r) • (V x a 



-V 



2m* (r) 



$A, g (r) = e x , q $A, ? (r) 



(8) 



where A runs over neutron and proton orbitals. 

For spherical nuclei, the single particle wave functions &\ t q can be fac- 
torized into a radial part and an angular part as (cf. Eq. (26) of Ref. 



0) 



where a is the spin and r the isospin, and the radial wave function satisfies 
the following equation: 



2m dr 2 2m r 2 9 



i>\,q( r ) = ^\,qi>\Ar) , (10) 



where (ft\ tq differs from the solution of Eq. (|10p . ip\ t q(r), by a normaliza- 
tion factor, 4>\ jq (r) = (m*(r)/m) 1 / 2 ^A,g( r )- The equivalent potential V q cq - is 
usually introduced for practical reasons, 



m*Jr) 



m 



V q (r) + U°°(r)(l ■ a) + 5 q , p V Co ui(r) 





- m*(rY 


+ 






m 



(11) 

where U q °(r) is the spin-orbit potential [23], Vcoul is the Coulomb potential, 
and 



V q {r) 



U, 



eff 



U^(r) + U co "(r) + U[ 
\ 2m* q {r) / h 2 
4 h 2 \2m*(r) 



12 



+ 



1 



h 2 



2 \2m*(r) 



+ 



H 2 



2m* (r] 



(12) 
(.13) 
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In the original work of Ma and Wambach [L5], the term U COTT (r) was de- 
rived directly from the Green's function with energy dependent self-energies 
while, in our approach, U corr {r) is derived from the new term @ in the 
EDF. A one-to-one correspondence between EDF and the Green's function 
approach is not possible. However, since the terms U co "(r) and (r) 
compensate each other in the Green's function approach [15], we want to 
reproduce the same behavior in the EDF. We obtain approximate compensa- 
tion between C7 corr (r) and U^ s (r) by imposing the following relation between 
the new coefficients: 

C p2(Vp)2 = 12 fm C T(Vp)2 . (14) 

We have investigated the sensitivity of the results to the value of the pro- 
portionality constant and we have checked that the reasonable values lie in 
the range 10-20 fm, after which the compensation is no longer efficient. 

The effects of the correction terms in @ will be analyzed in the next sec- 
tion. Notice that a surface-peaked effective mass could also be obtained from 
a modified Skyrme interaction. This different approach potentially leads to 
an improved agreement with experimental single particle energies |25j . Since 
the new term explored in Ref. [25] is simultaneously momentum and den- 
sity dependent, the functional obtained is quite different from Eq. (j4|): the 
number of terms is much larger and the correction to the effective mass is 
a polynomial in the density. It would be interesting to carry out a more 
detailed comparison of these two different approaches in a future study. 

2 Mean field properties 

In the following, we study the influence of the correction introduced by the 
new term Hq 011 for two representative nuclei: Ca and 208 Pb, using BSkl4 
interaction [26], which is adjusted to a large number of nuclei (2149). The 
effective mass in symmetric matter at saturation density is 0.8m and the 
isospin splitting of the effective mass in asymmetric matter qualitatively 
agrees with the expected behavior deduced from microscopic Brueckner- 
Hartree-Fock theory [27j . In the following, we study the effect of the correc- 
tion term on top of BSkl4 interaction without refitting the parameters. The 
refit is in principle necessary since the correction term impacts the masses 
and changes the single particle energies. In this first exploratory work, we 
improve the level density and discuss the effects on other quantities like the 
pairing, the entropy and the specific heat. The correction term (J4|) could 
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however potentially bring a better systematics in the comparison to experi- 
mental single particle centroids. This will be studied in a future work. 

The neutron and proton effective masses are plotted in Fig.[T]for different 
values of the coefficient Cg (Vp)2 = 0, -400, -800 MeV fm 10 . Notice that the 
case Cq^ p ^ = is that of the original Skyrme interaction BSkl4. Increasing 
|Co {Vp)2 | from to 800 MeV fm 10 , we observe an increase of the effective 

mass m* jvn at the surface, which produces a peak for large values of C T ^ P ^ , 
while at the center of the nucleus the effective masses get closer to the values 
obtained without the correction term. Fig. [T]can be compared with Fig. 1 of 
Ref. [15]. Due to the different surface-peaked functions (here (Vp) 2 instead 
of a single V dependence in Ref. [E]), the width of the effective mass at 
the surface is larger in Ref. |15j than in the present work. For values of the 



coefficient \C T ^ P ^ | larger than 800 MeV fm 10 , the potentials U COTT (r) and 
Uq S (r) (Eq. ([7D, (f!3])) induce large gradients of the mean field (JT2J) in a tiny 
region close to the surface of the nuclei, which in turn produce an instability 
in the HF iterations. 

A surface-peaked effective mass has also been deduced from the particle- 
vibration coupling within the HF+RPA framework [20j. In such an ap- 
proach, the effective mass is shown to be peaked not only at the surface 
of the nuclei, but also in a window around the Fermi energy of ± 5 MeV. 
Such an energy dependence could not be implemented straightforwardly in 
the EDF framework. We could however evaluate the state-averaged effective 
mass, (m*/m}\, defined as: 

KMa = fdr ^ 9 ( r )^^ A)? (r) , (15) 

where the index A stands for the considered state. The state-averaged effec- 
tive masses {m*/m)\ are represented in Fig. [2] as a function of the energy 

of the bound states, and for different values of the parameter Cq^ V ^ . As 

the value of the coefficient |Cq^ V/ ^ | gets larger, the state-averaged effective 
masses (m*/m)\ approach 1 around the Fermi energy while it get closer to 
the original effective mass for deeply bound states. Notice however quan- 
titative differences between Ca and Pb for states around the Fermi energy 
as well as deeply bound states. In conclusion, even if we did not introduce 
an explicit energy dependence of the surface-peaked effective mass, we still 
find that the expected behavior [20] of (m*/m) \ ELS £1 function of energy is 
qualitatively reproduced. 

In order to evaluate the impact of the new term on the density of states, 
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we display in Fig. [3] the neutron and proton number of states as a function 
of the excitation energy, defined as: 

N(E) = / dE' g(E') , (16) 
J o 

where g(E) is the density of states, 

g(E) = ^p-= Y, (2jx 2 + l)S(E-(e X2 -e Xl )) , (17) 

Ai<F,A 2 >F 

and e\ 1 (e\ 2 ) represent the single-particle energies below (above) the Fermi 
surface. The expected relation between the surface-peaked effective mass 
and the density of states is clearly shown in Fig. [3j the number of states 
at given excitation energy increases as the coefficient Cq^ p ^ goes from 

to -800, meaning that the density of states also increases as |Cq^ Vp ^ | gets 
larger. 

Let us now discuss qualitatively the impact of the correction term on 
global properties of nuclei starting with the density profiles. The neutron 
and proton densities are shown in Fig. |4l Since for 40 Ca the neutron density 
is very similar to the proton one, we have chosen to represent only the neu- 
tron one. Since the number of particles has to be conserved, lower values of 
the density in the bulk of nuclei for larger values of \Cj^ p ^ | are compen- 
sated by a slight increase of the size of the nucleus. These small differences 
in the density profile influence the charge root mean square radius r c h (see 

Table [T]), which slightly increases as the value of the parameter \Cq^ p ^ \ 
gets larger. Moreover, because of the isoscalar nature of the correction (j3J, 
neutrons and protons are affected in an identical way, as one can see from 
the constant value of the neutron skin radius, r s ki n , given in Table HJ 

Let us now analyze the influence of Hq ovv at the level of the mean field 
V q (r) defined in Eq. (I12p . The different components of the central part of 
the mean field, Uq ky (r), U^ s (r), and U corr (r), see Eqs. ([7]) and (fl~3j) . are 
represented in Fig. [5] for the neutrons, the protons, and for 40 Ca and 208 Pb 
nuclei. The value of the coefficient is fixed to be Cq^ p ^ = —400 MeV 
fm 10 . As expected, there is a reasonable compensation between Ug S (r) and 
U corr (r). As a consequence, the mean field V q {r) is nearly not affected by the 
presence of a surface-peaked effective mass for the values of the coefficient 
£^(Vp) c hosen in the domain going from to -800 MeV fm 10 , except close 
to the surface where a small change in the slope is observed. We have then 
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shown that in the EDF framework, the correction term (J3|) reproduces the 
result obtained in Ref. [15] . 

We now compare our results to the recent ones of Zalewski et al. |28y29j 
where correction terms such as as well as others have been studied. The 
main differences between our approach and that of Refs. [2"5]I59] are: (i) the 
moderating term is not present in Refs. [28,29], (ii) the effective mass in the 
bulk of nuclei is close to 0.8m in our case, while in the Refs. [28J29J it assumes 
the value 1 since they started with SkX-Skyrme interaction [31J , (Hi) we 
have not refitted the parameters of the interaction contrarily to Refs. [28} 
[29] ; indeed, at variance with our approach, the functional is readjusted in 
Refs. [28,29j such that the condition 

[ dr po^nm =1 (18) 

J A m v ' 

is satisfied. An important dependence of the spin-orbit splitting and of the 

centroids as a function of the coefficient of the correction term have been 

observed in Ref. [28]. In our case, as shown in Figs. [6] and [7] for both Ca 

and 208 Pb nuclei, we do not find such an important effect. Only a weak 

p 2 (Vo) 2 

dependence on the coefficient Cq of the spin-orbit splitting is observed 
and an almost independence of the spin-orbit centroids. This contradiction 
between the two results might come from the readjustement of the Skyrme 
parameters (point (Hi) mentioned above) performed in Ref. [29]. In the 
latter article, the spin-orbit interaction is not changed, but the parameters 
of the Skyrme interaction are changed, which might induce a modification 
of the density profile and, therefore, of the spin-orbit splitting. In fact, in 
Ref. [29] another procedure is adopted, the non-perturbative one, and the 
behavior of the the spin-orbit splitting and of the spin-orbit centroids is quite 
different from that shown in Ref. [28J. The non-perturbative prescription of 
Ref. [29] shows almost no change of the spin-orbit splittings and centroids 
with respect to the strength of the surface-peaked effective mass. This shows 
that the spin-orbit splittings and their centroids might not be impacted by 
the presence of a surface-peaked effective mass in a direct way, but eventu- 
ally, indirectly, through the readjustment procedure of the functional. 

Finally, we have studied how the binding energy varies as a function of 
the correction term (|4|). The results are presented in Table [2j The binding 

energy of 40 Ca and 208 Pb increases as the parameter \Cj^ p ^ | increases. We 
can therefore expect that the readjustment of the parameters of the Skyrme 
interaction shall essentially make the interaction more attractive. 
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3 Pairing properties 



Most of the nuclei are superfluid and it could be shown that in the weak cou- 
pling limit of the BCS approximation, the pairing gap at the Fermi surface 
Ap and the pairing interaction v pa i r are related in uniform matter through 
the relation [30J: 

Ap ~ 2ep enp[2/ (NoVpair)], (19) 

where ep is the Fermi energy, iVo = m*kp/(h 2 ir 2 ) is the density of states at 
the Fermi surface. Then a small change of the effective mass m* can result 
in a substantial change of the pairing gap. 

In the following, we will consider 120 Sn because it is an excellent candi- 
date to study pairing correlations [23J: 120 Sn is spherical and only neutrons 
are participating to the S-wave Cooper-pairs. An accurate description of 
the pairing properties can be obtained already at the level of the spherical 
HF+BCS framework |23| . In this section, we study qualitatively the relation 
between the increase of the effective mass at the surface and its consequences 
on the pairing properties, both at zero and at finite temperature. 

We adopt a density-dependent contact interaction v nn given by [33] : 

(k\v nn \k') =v g( P )9(k,k>) , (20) 

where the strength vo of the pairing interaction is adjusted to obtain an av- 
erage pairing gap equals to 1.3 MeV in 120 Sn. The factor g(p) in Eq. ([20]) is a 
density-dependent function (see below) and 8(k,k') is the cutoff introduced 
to regularize the ultraviolet divergence in the gap equation. We choose for 
9(k,k') the following prescription: 9(k,k') = 1 if E^^E^i < E c , otherwise it 

is smoothed out with the Gaussian function exp ^— [(Ex — E c )/a] 2 ^j. Here- 
after, we choose E c = 8 MeV and a = 1 MeV. Notice that to be compatible 
with HFB calculations, the cutoff is implemented on the quasiparticle en- 
ergy, E\ s = [(e\ t q— /ig) 2 +A 2 J 1 / 2 , where e\ is the HF energy, \i the chemical 
potential, and A\ A the average pairing gap for the state A (see Eq. ([26j) ) . 
The density-dependent term g(p) is simply defined as: 

g(p) =l-v~, (21) 
Po 

where ij designates the volume (rj = 0), mixed (rj = 0.5), or surface (r) = 1) 
character of the interaction and po = 0.16 fm -3 is the saturation density of 
symmetric nuclear matter. The isoscalar particle density, p = p n + p p , is 
defined as: 

P^ r ) = h + !) KiQ- ~ /A«) + U U^] l^»! 2 , (22) 

A 
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where the u\ A and the v\ >q are variational parameters. The represent the 
probabilities that a pairing state is occupied in a state (A, q), = 1 — v\ , 
and fx ;Q is the Fermi function for the quasiparticle energy E\ tq [53] : 

fx > q = i + e E Kq /k B T ' ( 23 ) 

where fee is the Boltzmann constant. 
The local pairing field is then given by: 



A 9 (r) = a ^(p)p fl (r) , (24) 



where /5 g is the abnormal density defined as: 

P 9 (r) = J>j A + l)«A, g «A, g (l " 2/A, 9 )|<AA l(? (r)| 2 . (25) 



47T 

A 



The Eq. (|24j) is the self-consistent gap equation that should be solved con- 
sistently with the particle conservation equation (|22p . The average pairing 
gap for the state A used in the definition of the cutoff is defined as: 

A Aj(? = j dr|</> A »| 2 A g (r). (26) 



In the HF+BCS framework, Eqs. (|22p and (|24p are solved at each iter- 
ation. The number of iterations performed depends on the convergence of 
the pairing gap equation (f24"l) : it goes from about 200 up to about 1000 near 
the critical temperature. 



3.1 Pairing properties at T = 

In the following, we study the influence of the correction term Q on the 
pairing properties at T = 0. A realistic calculation for nuclei shall treat 
consistently the pairing interaction in the particle-particle channel and that 
in the particle-hole channel. The bare interaction in the particle-particle 
channel shall then be replaced by the induced one which accounts for a 
50% correction [35] . It is however not our intention in the present work 
to investigate this question. We want, at a simpler level, to clarify the 
role of the correction term Q on the pairing properties, and we will show 
that there is indeed a correlation in space between the enhancement of the 
effective mass and that of the probability distribution of the Cooper pairs. 
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In order to reproduce the value of the average gap A n = 1.3 MeV in 
120 Sn, we have adjusted Eq. (|24[) for three kinds of pairing interactions (vol- 
ume, mixed, and surface), for the functional without the correction term @. 
We obtain for the values (r),v ): (0; -259 MeV fm 3 ), (0.5; -391 MeV fm 3 ), 
(1; -800 MeV fm 3 ). In Table Owe report the values for the average neutron 
pairing gap A n , calculated as the average gap over the abnormal density, 

A n = i / drp n (r)A n (r) , (27) 



N 

where N = fdr p n (r), and for several values of the coefficient Cjf^ p ^ . In 
Table [3l it is shown that the effect of the correction term @ on the pairing 
gap is non-negligible. The pairing gap is increased by 200 to 700 keV as the 
coefficient \Cq^ p ^ | gets larger. An important dependence with respect to 
the kind of the pairing interaction (volume, mixed, surface) is also observed. 
The largest effect of the correction term Q is obtained for the surface pairing 
gap. From the results presented in Table [3] we notice that the correction 
term @ in the functional has an important influence on the average pairing 
gap. 

In order to study the effect of the correction term on the pairing prop- 
erties in a more realistic case, we have chosen to renormalize the pairing 
interaction, for each value of the coefficient Cq^^ , in such a way to get 
always at zero temperature A n = 1.3 MeV. We obtain, at T = 0, for 
c r(v P ) 2 = _ 40Q MeV fm i 0; the valueg ( V}VQ y m ( 0; .248 MeV fm 3 ), (0.5; 

-362 MeV fm 3 ), (1; -670 MeV fm 3 ), and, for = -800 MeV fm 10 , the 

values (v,v ): (0; -235 MeV fm 3 ), (0.5; -337 MeV fm 3 ), (1; -593 MeV fm 3 ). 
We represent in Fig. [8] the pairing field and the probability distribution of 
Cooper-pairs defined as 

p(r) = -4irr 2 p(r) , (28) 

versus the radial coordinate for different values of the parameter Cq^ p ^ . 
On the left panels, we display the pairing field profiles, which do not change 
very much with the coefficient Cj^ p ^ ; this is due to our renormalization 
procedure, which requires the average pairing gap to be 1.3 MeV. The right 
panels of Fig. [8] clearly show the correlation in space between the enhance- 
ment of the effective mass and that of the probability distribution of the 
Cooper pairs, i.e. the enhancement of the probability distribution is located 
where the effective mass is surface-peaked. 

In conclusion, it has been shown in this section that the effect of the 
surface-peaked effective mass on the pairing properties is non-negligible. 
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As for approaches where the pairing interaction is empirically adjusted on 
some nuclei, we absorbed this effect through a renormalization of the pairing 
interaction. For non empirical approaches where the pairing interaction 
is not adjusted on nuclei properties but directly to that of the bare 1 So 
potential [36H50] , the enhancement of the average pairing gap induced by the 
surface-peaked effective mass shall be treated consistently with the induced 
pairing interaction [35j . 

3.2 Pairing properties at finite temperature 

The density of states around the Fermi energy shall influence the temperature- 
related quantities such as the entropy and the specific heat [5], and at the 
same time the density of states is also affected by the temperature [6| . l42 |. l43j . 
We explore the effect of the new term on 120 Sn in the framework of HF+BCS 
at finite temperature, using the pairing interaction where the strength has 
been renormalized at T = for each value of the coefficient Cq^ V ^ . In Fig. [9] 
we show the neutron pairing gap as a function of temperature, for the three 
different kinds of pairing interactions (volume, mixed, surface). The results 

on the left are obtained keeping the coefficient Cq^ V ^ constant, while, on 

the right, we plot the results obtained letting Cq^^ vary exponentially 
with temperature, according to the following relation: 

^(Vp) 2 (r) = C r(V^ xe -T/T ; (29) 

with To = 2 MeV. The choice of this kind of temperature dependence relies 
on the work by Donati et al. [6J, where a study of the w-mass in the frame- 
work of QRPA for temperatures up to some MeV was carried out on some 
neutron rich nuclei. The variation of m u with respect to temperature was 
parameterized with an exponential profile and the typical scale of the varia- 
tion of m u with temperature was found to be around 2 MeV. Notice that in 
semi-infinite nuclear matter, the typical scale was found to be 1.1 MeV [43] . 
The reduction of the scale might be due to the semi-infinite model which is 
still far from a realistic finite nucleus case. 

We observe in all cases in Fig. [9] the typical behavior associated to the 
existence of a critical temperature T c after which pairing correlations are 
destroyed In particular, as expected, T c is not modified by the new 

term for the case of constant coefficient, since we absorbed the effect of the 
new term through the renormalization procedure. Instead, the effect of a 
temperature dependent coefficient Cq^ V ^ (T) is to reduce the critical tem- 
perature; indeed, the chosen dependence (|29l) shifts the critical temperature 
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of ~ 40 keV in the case of volume interaction and of ~ 60 keV in the case 
of surface interaction. The relation T c ~ A n (T = 0)/2 is still verified; more 
precisely, the ratio T c /A n (T = 0) is ~ 0.55 (for volume interaction) and 

~ 0.57 (for surface interaction) for the case of T-independent Cq^ p ^ , while 
with the T-dependent prescription (|29p it varies from ~ 0.55 (for volume 
interaction) to ~ 0.51 (for surface interaction). In conclusion, we remark 
that through its temperature dependence, the surface-peaked effective mass 
has an effect on the critical temperature, that could also be extracted exper- 
imentally |44y45j . This perspective motivates the application of the present 
work to realistic cases at finite temperature. 

Finally, in Fig. [10] we show the total entropy and specific heat as a 
function of temperature for volume and surface interaction, and for three 
values of Cq {Vp)2 (0, -800 MeV fm 10 , and Cq [ ^ p)2 (T) for the case -800 MeV 
fm 10 ). The entropy of the system is calculated as: Stot = S n + S p , being: 

S q = -k B \f*a ln (^<?) + ( X - ln (! " A,,)] , (30) 

A 

where f\q is defined as in Eq. ([23]) . The specific heat is then defined as: 

Cv = T^. (31) 

We observe the change of the slope in the entropy curve, which causes the 
discontinuity in the specific heat in correspondence of the critical tempera- 
ture. In agreement with the previous results, T c is shown to be modified by 
a temperature-dependent Cq^ p ^ , the effect is stronger in the case of surface 
interaction, and the temperature dependent coefficient acts as to reduce T c . 



4 Conclusions 

In this paper we have studied the influence of the correction term Q on 
various nuclear properties. The isoscalar correction term Q has been shown 
to produce a surface-peaked effective mass in the nuclei under study ( 40 Ca 
and 208 Pb), without modifying significantly the mean field profiles. The 
increase of m* jm at the surface is up to about 1.2 - 1.3 for the maximum 
value of the strength of the correction we used. As the effective mass gets 
enhanced at the surface of nuclei, the density of states increases. Then, we 
have studied the impact of such a term on the neutron pairing gap in the 
semi-magic nucleus 120 Sn, within an HF+BCS framework, and it turned out 
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that its effect is non-negligible; if the pairing interaction is not renormalized 
consistently with the new term, the average gap increases from 200 keV to 
700 keV under variation of the strength of the correction term and depending 
on the volume/surface character of the interaction. In a recent work [36], it 
has been stressed that the surface enhancement of the pairing field induced 
by the PVC might also play a role on the size of the Cooper-pairs at the 
surface of nuclei. In uniform matter, the coherence length is indeed inversely 
proportional to the pairing gap. The surface enhancement of the pairing field 
could then make the Cooper-pairs smaller at the surface of nuclei. It would 
be interesting to investigate whether this interesting feature of the pairing 
correlation might be probed experimentally, for example by pair transfer 
reaction mechanism. Finally, we have explored some finite temperature 
properties in 120 Sn, within a HF+BCS framework at finite temperature. 
We observed for the neutron pairing gap that the critical temperature at 
which pairing correlations vanish is shifted if a T-dependence in the new 
coefficient is considered. As a consequence, the entropy and specific heat 
profiles are affected by the introduction of the new term. 

In the future, we shall go on with a global refitting of all the parameters 
of the functional, including the new correction term @. It would be very 
instructive to know whether both the masses and the single particle levels 
could be improved in such a way. 

A. F. F. and J. M. would like to thank E. Khan and N. Sandulescu for 
their help in checking our HF+BCS results, as well as the fruitful discussions 
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Table 1: Charge rms radius, r c h, and neutron skin radius, r s k in , for °Ca 
and 208 Pb and for different values of the coefficient . The charge rms 



radius is calculated according 


; to Eq. (110) in Ref. [23 . 
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r c h 
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[MeVfm 10 ] 
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0.16 
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-600 
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0.16 



Table 2: Binding energy per nucleon in 40 Ca and 208 Pb (in MeV) for different 



values of the coefficient Cq^ p ^ . 




40 Ca 


208p b 
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-200 
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-7.983 


-400 
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-600 
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-7.849 


-800 


-8.227 


-7.801 
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Table 3: Average neutron pairing gap A ra for 120 Sn for different values of 
the coefficient CjS^ p ^ . 
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Figure 1: m*/m as a function of radial coordinate for °Ca and 208 Pb, for 



c; 



r(Vp) 2 



0, -400, -800 MeV fm 10 . 
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Figure 2: Energy dependence of the expectation value of (m*/m)\ for bound 

states for 40 Ca and 208 Pb, for Cq (Vp)2 = 0, -400, -800 MeV fm 10 . The arrows 

indicate the position of the Fermi energy for Cq {Vp)2 = MeV fm 10 defined 
as the energy of the last occupied state. The correction induced by the 
surface-peaked effective mass produces an increase of the Fermi energy by 
400 keV at most. 
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Figure 4: Neutron density as a function of radial coordinate for 40 Ca and 
neutron and proton densities as a function of the radial coordinate for 208 Pb, 
for Cg (Vp)2 = 0, -400, -800 MeV fm 10 . 
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Figure 5: Central part of the neutron and proton mean field as a function 
of radial coordinate for 40 Ca and 208 Pb, for = -400 MeV fm 10 . 
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Figure 6: Spin-orbit splitting and centroids for 40 Ca as a function of the 
The experimental values are taken from |32j. 



coefficient Cq^ Vp ^ 



24 



neutrons 



protons 



5 7 
a 3 

■I 4 



> 

-a 

'8 

c 

u 
u 

o 



.2g 

2f 



3p 7" 



•■•3d- 

~3g" 
2i 

A 

3 £. 



2f 



200 



lh 























lg 


















"5 










A 


2d 

II 






















2f 








■ - 


■ 


lh 

■i----- 


.-=.:K--- 


-.:-=.:A-- 




._.?d.. A 



lg 



400 

, t(gradp) 



600 



800 



200 



7 
6 
5 
4 
3 
2 
1 



-2 

-4 

-6 

-8 

-10 

-12 

-14 



400 

i ^(gradp)" 



600 



800 



Figure 7: Spin-orbit splitting and centroids for 208 Pb as a function of the 



coefficient Cq^ Vp ^ 



The experimental values are taken from [24J. 
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Figure 8: Neutron pairing gap (on the left) and probability distribution of 
the Cooper pairs (on the right) for 120 Sn as a function of radial coordinate, 

for Cq (Vp) = 0, -400, -800 MeV fm 10 , and for different types of pairing 
interaction (volume, mixed, surface). 
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Figure 9: Neutron pairing gap for 120 Sn as a function of temperature for 
different values of the coefficient Cq^ V ^ . The results on the left are ob- 
tained keeping Cq (Vp)2 constant {C^^ = 0, -400, -800 MeV fm 10 ), while 

the results on the right are obtained using for Cj^ p ^ the prescription in 
Eq. (EH]). 
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Figure 10: Total entropy Stat and specific heat Cy (in units of the Boltzmann 

constant), for 120 Sn as a function of temperature, for C^ p ^ = 0, -800 MeV 
fm 10 , -800 e~ T l T ° MeV fm 10 , and for volume (on the left) and surface (on 
the right) pairing interaction. 
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